<html><head><title>Comparing centiles from different GAMLSS models</title>
<meta http-equiv="Content-Type" content="text/html; charset=">
<link rel="stylesheet" type="text/css" href="Rchm.css">
</head>
<body>

<table width="100%"><tr><td>centiles.com(gamlss)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
<param name="keyword" value="R:   centiles.com">
<param name="keyword" value=" Comparing centiles from different GAMLSS models">
</object>


<h2>Comparing centiles from different GAMLSS models</h2>


<h3>Description</h3>

<p>
This function compares centiles curves for more than one GAMLSS objects. 
It is based on the <code>centiles</code> function. 
The function also tabulates the sample percentages below each centile curve (for comparison with the model percentages given by the argument <code>cent</code>.)
A restriction of the function is that it applies to models with one explanatory variable only.
</p>


<h3>Usage</h3>

<pre>
centiles.com(obj, ..., xvar = NULL, cent = c(0.4, 10, 50, 90, 99.6), 
             legend = TRUE, ylab = "y", xlab = "x", xleg = min(xvar), 
             yleg = max(obj$y), xlim = range(xvar), ylim = NULL, 
             no.data = FALSE, color = TRUE, plot = TRUE)
</pre>


<h3>Arguments</h3>

<table summary="R argblock">
<tr valign="top"><td><code>obj</code></td>
<td>
a fitted gamlss object from fitting a gamlss continuous distribution </td></tr>
<tr valign="top"><td><code>...</code></td>
<td>
optionally more fitted GAMLSS model objects</td></tr>
<tr valign="top"><td><code>xvar</code></td>
<td>
the unique explanatory variable </td></tr>
<tr valign="top"><td><code>cent</code></td>
<td>
a vector with elements the % centile values for which the centile curves have to be evaluated</td></tr>
<tr valign="top"><td><code>legend</code></td>
<td>
whether a legend is required in the plot or not, the default is <code>legent=TRUE</code> </td></tr>
<tr valign="top"><td><code>ylab</code></td>
<td>
the y-variable label </td></tr>
<tr valign="top"><td><code>xlab</code></td>
<td>
the x-variable label </td></tr>
<tr valign="top"><td><code>xleg</code></td>
<td>
position of the legend in the x-axis</td></tr>
<tr valign="top"><td><code>yleg</code></td>
<td>
position of the legend in the y-axis</td></tr>
<tr valign="top"><td><code>xlim</code></td>
<td>
the limits of the x-axis</td></tr>
<tr valign="top"><td><code>ylim</code></td>
<td>
the limits of the y-axis</td></tr>
<tr valign="top"><td><code>no.data</code></td>
<td>
whether the data should plotted, default <code>no.data=FALSE</code> or not <code>no.data=TRUE</code> </td></tr>
<tr valign="top"><td><code>color</code></td>
<td>
whether the fitted centiles are shown in colour, <code>color=TRUE</code> (the default) or not <code>color=FALSE</code> </td></tr>
<tr valign="top"><td><code>plot</code></td>
<td>
whether to plot the centiles</td></tr>
</table>

<h3>Details</h3>




<h3>Value</h3>

<p>
Centile plots are produced for the different fitted models and the sample centiles below each centile curve are printed</p>

<h3>Warning</h3>

<p>
This function is appropriate only when one continuous explanatory variable is fitted in the model
</p>


<h3>Note</h3>




<h3>Author(s)</h3>

<p>
Mikis Stasinopoulos <a href="mailto:d.stasinopoulos@londonmet.ac.uk">d.stasinopoulos@londonmet.ac.uk</a> and Bob Rigby <a href="mailto:r.rigby@londonmet.ac.uk">r.rigby@londonmet.ac.uk</a>
</p>


<h3>References</h3>

<p>
Rigby, R. A. and  Stasinopoulos D. M. (2005). Generalized additive models for location, scale and shape,(with discussion), 
<EM>Appl. Statist.</EM>, <B>54</B>, part 3, pp 507-554.
</p>
<p>
Stasinopoulos D. M., Rigby R.A. and Akantziliotou C. (2006) Instructions on how to use the GAMLSS package in R.
Accompanying documentation in the current GAMLSS  help files, (see also  <a href="http://www.gamlss.com/">http://www.gamlss.com/</a>).
</p>


<h3>See Also</h3>

<p>
<code><a href="gamlss.html">gamlss</a></code>, <code><a href="centiles.html">centiles</a></code> , <code><a href="centiles.split.html">centiles.split</a></code>
</p>


<h3>Examples</h3>

<pre>
data(abdom)
h1&lt;-gamlss(y~cs(x,df=3), sigma.formula=~cs(x,1), family=BCT, data=abdom)
h2&lt;-gamlss(y~lo(x,span=0.4), sigma.formula=~lo(x,span=0.4), family=BCT, data=abdom ) 
centiles.com(h1,h2,xvar=abdom$x)
rm(h1,h2)
</pre>



<hr><div align="center">[Package <em>gamlss</em> version 1.6-5 <a href="00Index.html">Index]</a></div>

</body></html>
